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THE PRESSURE FUNCTION FOR PRODUCTS 
OF NON-NEGATIVE MATRICES 

DE-JUN FENG AND KA-SING LAU 



Abstract. Let (S^, a) be a subshift of finite type and let M(x) be a continuous func- 
tion on Y^a taking values in the set of non- negative matrices. We extend the classical 
scalar pressure function to this new setting and prove the existence of the Gibbs mea- 
£/3 ' sure and the differentiability of the pressure function. We are especially interested on 

the case where M(x) takes finite values Mi, • • • ,M m . The pressure function reduces to 
P(q) '■= linin^oo — log ^2 jgyj A ll-^/ll 9 - The expression is important when we consider 
the multifractal formalism for certain iterated function systems with overlaps. 

1. Introduction 

> 

Let a be the shift map on S = {1, 2, • ■ • , m} N , m > 2. As usual £ is endowed 
with the metric d(x, y) = m~ n where x = (xf.), y = {yk) and n is the smallest of the 
k such that Xk ^ yk- Given an m x m matrix A with entries or 1, we consider the 
subshift of finite type (S^, a) (see 0). We shall always assume that A is primitive. 

Suppose M is a continuous function on T,a taking values in the set of all non- 
negative d x d matrices. For q G R, we define the pressure function P(q) of M 
by 

P(q) = lim - log V sup \\M(x)M(ax) . . . Mta^x)^, (1.1) 

n— *rsn n f J _r t t 



where n denotes the set of all admissible indices of length n over {1, . . . , m}; for 
J — ji ' " ' jn G Sa,h, [J] denotes the cylinder set {x = (xi) G XU : X{ = ji, 1 < 
i — || ' || denotes the matrix norm denned by ||S|| := 1*51, 1* = (1, 1, ... , 1). 
By using a subadditive argument, it is easy to show that for q > 0, the limit in 
the above definition exists. With some additional conditions on the matrices (see 
Theorem 1.1), the limit exists for q G R. 



Key words and phrases: Pressure, Product of matrices, Gibbs measures, Iterated function sys- 
tems, Hausdorff dimension, Multifractal formalism. 
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The pressure function of the scalar case (i.e., M(x) = where <f>(x) is a real 



valued function called the potential of the subshift) has been studied in great detail 
in statistical mechanics and dynamical systems in conjunction with the Gibbs mea- 
sure, the entropy and the variational principle (c.f., e.g., [B], [P], [R]); it has also 
been used to study the multifractal structure of the self-similar (or self-conformal) 
measures generated by iterated function systems (IFS) with no overlap (the open 
set condition) ([MU], [FL]). By identifying with the symbolic space, such self-similar 
measure \x is actually a Gibbs measure and the pressure function is directly related 
to the scaling spectrum of \i [FL, Theorem 3.3]. In all the above cases, the pressure 
functions under consideration are differentiable (actually real analytic). This prop- 
erty is essential to investigate the phase transition in thermodynamics and for the 
multifractal formalism in the dimension theory of fractals. 

In the recent investigation of the self-similar measures generated by iterated func- 
tion systems with overlaps, it is seen that in many interesting cases, such measure 
fi can be put into a vector form with a new non-overlapping IFS and with matrix 
weights ([LN1,2], [LNR], [Fe], [FeO]). In this way the validity of the multifractal for- 
malism depends on the differentiability of the pressure function P(q) in (1.1) (more 
precisely (1.4) in the following) [LN2]. In another direction, the expression of the 
matrix product in (1.1) also appears in the study of the scaling functions in wavelet 
theory (the matrices are allowed to have negative entries) in the form of L 9 -joint 
spectral radius and the L 9 -Lipschitz exponent ([DL1, 2], [LM]); the problem of dif- 
ferentiability of the P(q) also appears there. So far there is no general theorem to 
guarantee this fact other than some special cases (e.g., [LN1], [FLN], [Fe], [FeO], 



The main purpose of this paper is to consider the pressure functions and the Gibbs 
measures for the products of matrices. We first study the case that the matrices 
M(x), x G XU are positive, we prove the following fundamental theorems. 

Theorem 1.1. Suppose that M is a Holder continuous function on T, A taking values 
in the set of positive d x d matrices. Then for any q G K, there is a unique a- 
invariant, ergodic probability measure ji q on of which one can find constants 
C 1 > 0, C 2 > such that 



[DL2]). 



Ci< 




< C 2 



(1.2) 



exp(-raP(g)) • \\M(x)M(ax) . . . M(a n ~ 1 x)\\ q 



for any n > 0, J 



G T, An and x G [J]. 
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The above measure fi q is called the Gibbs measure associated with M and q. 
We remark that the theorem generalizes the classical existence result of the Gibbs 
measure for a real- valued M(x) (see [B, §1.4]). The positivity of the matrices is used 
to yield the follow simple estimate (Lemma 2.1) 

\\M(x) ■ ••M(o n+t ~ 1 x)\\ « \\M(x)---M(a n - 1 x)\\ ■ \\M(a n x) ■ ■ ■ M(a n+i ~ 1 x)\\. (1.3) 

By using this we can apply a technique of Brown, Michon and J. Peyriere [BMP] 
and Carleson [C] to construct a certain ergodic measure which is the Gibbs measure 
fi q . The fig has the following quasi- Bernoulli property (Heurteaux [H]): there exists 
C > such that for any n, k G N with I G £U, n , J e ^U,* and IJ G IU,n+fc 

C-V 9 ([/])^([J]) < ^([IJ]) < C/i 9 ([/])// ? ([J]) (1.4) 
This together with a result of Heurteaux [H] imply 

Theorem 1.2. Under the condition of Theorem \l.\ P(q) is differentiable for q ^ 0. 
As an application, we let 

log||M(s)M(<rs)---M(<r"- 1 (s)|| 
-Mm := iG La : lim = a > . 

k n— >oo n 

We prove the following dimension formula 

Theorem 1.3. Under the same assumption of Theorem [771| ; we /love /or any a = 
P'(q), q^o, 

dim H E(a) = (—ocq + P(q)) 

iogm 

where dimn denotes the Hausdorff dimension. 

The above theorems depend very much on the positivity of the matrix-valued 
M(x). In order to extend them to nonnegative matrix- valued functions, we have to 
impose more conditions on M(x): 

(HI) M(x) = Mi if x G [i], i = 1, • • • , m\ 

(H2) M is irreducible in the following sense: there exists r > such that 
for any i,j G {1,2, . . . ,m}, 

r 

where T>A,k;i,j denotes the set of all K G XU,fc such that iKj G TiA,k+2- 
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We see that under the assumption (HI), the pressure function in (1.1) can be re- 
written as 

P(q) = lim - log V \\Mj\\ q . (1.6) 



A,n 



where Mj = Mj. ■ ■ -Mj n . If = £ is the symbolic space with a full shift, then 
(H2) is equivalent to ELi( M i + " " " + M ™) k > 

In this new setting, we use (H2) to adjust (1.3) and the required lemmas, the 
Gibbs measure \i q is shown to exist for q > 0. This time \x q only satisfies fi q ([IJ]) < 
C n q ([I}) [i q ([J]) instead of (1.4); nevertheless we can still prove the differentiability of 
P(q), q > as in Theorem 1.2. Theorem 1.3 can be adjusted likewise (see Theorem 
3.3 and Theorem 3.4). As an application, the first author proves the smoothness 
of the L 9 -spectrum (q > 0) and the multifractal formalism for a class of self-similar 
measures with overlaps (including the Bernoulli convolutions associated with Pisot 
numbers) in a forthcoming paper [Fe2]. 

For the organization of the paper, we prove the above results for the positive 
matrix-valued functions in Section 2. In Section 3, we modify the proofs for the 
non-negative matrix- valued functions with (HI) and (H2). In Section 4, we give an 
illustration of reducing an IFS with overlap to a vector-valued IFS with no overlap, 
and the pressure function in (1.6) arises. We also give some remark on the theorems 
and raise a few unsettled problems. 
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it to the present form which can be appealed to more general situation. They also 
thank E. Olivier for introducing the multifractal results of [[BMP], |H| and O. Stenflo 



for reading the manuscript carefully and suggesting some improvements. 



2. Positive Matrices 

In this section we assume that M is a Holder continuous function on taking 
values in the set of all positive d x d matrices. 

For any two families of positive numbers {ai} ie j, {bi} ieI , we write, for brevity, 

w bi to mean the existence of a constant C > such that C _1 aj < bi < Ca>i for 
all i a I\ a,i ^ bi to mean the existence of a constant C > such that < Cbi for 
all i G X and a« ^ bi means 6j =^ a^. 

We start with a simple lemma: 
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Lemma 2.1. For any x E Ha, n, £ E N, 

\\M(x)---M(a n+e - 1 x)\\ « HM^-.-M^- 1 ^!! • ||M((T n x) • • • M((j" + ^ _1 x) || 

(the involved constant in zs independent of n,£ and x). 

Proof. It is clear that 

\\M{x) ■ ■ ■ M(or n+l ~ l x) || < \\M{x) ■ ■ ■ M(a n - l x) \\ ■ \\M(a n x) ■ ■ ■ M^+^x) \\ . 

To prove the reverse inequality, we observe that M is positive and continuous, there 
is a constant C > such that 

maxj j Mi j [x) 

This implies that M(x) > ^-EM(x) (A > B means that A it j > B it j for each index 
(hj)) where E = {Eij)i<ij<d is the matrix whose entries are all equal to 1. Let 1 
be the ci-dimensional column vector each coordinate of which is 1. Then 

\\M(x) ■ ■ ■ M(a n+e ~ l x) || > \\M(x) ■ ■ ■ M(a n - l x)^-EM(a n x) ■ ■ ■ M{a n+i ~ l x) \\ 

= %\\M(x) ■ ■ ■ M((r n_1 x)lTM(fT n a;) • • • M(a n+ ^ l x)\\ 
a 

= ^-\\M(x)---M(a n ' 1 x)\\ ■ \\M(a n x)---M(a n+e ~ 1 x)\\. 

□ 



We define 



s n (I,q) = sup \\M(x)M(ax) . . . M(a n ~ 1 x)\\ q V/ E E A , 
xe[i] 



and 



Lemma 2.2. For a fixed q E R ; 

s n (I, q) » \\M(x)M(ax) . . . M^x)^ V / E E A , n , x E [/]. 
Proof. For any fiGN, define 

Vn = sup = I e £a,„, ije [/], 1 < i, j < d j . (2.2) 



Since each Mij is positive and Holder continuous, we have |log?y n | < Cm an for 
some C > and < a < 1. It follows easily that r] := n^Li Vn < 00 an d hence for 
x G [/], 

\\M{x)...M{a n - 1 x)\\ q < s n (I,q) < rj^ \\M(x) . . . M(a n ^ 1 x) \\ 9 . 

□ 

We have assumed that A is primitive, there is an integer p > such that A p > 0. 
This implies that for any I G S^^, J G there exists K G E^ iP such that 

IK J G S^ >n+ £ +p . 

Lemma 2.3. Let p be such that A p > 0. Then for a fixed q G M, 
(i) s/(g) w s*_ p (g) /or a// £ > p; 

(ii) For / G Ea,„, ? > P, J2j s n+t(U, q) ~ s n+e(JI, <?) ~ s n (I, g)s*(g) 
where the first (second) sum is taken over all J G E^ such that I J G Y,A,n+e 
(JI E T,A, n +e respectively ); 

( iu ) iKj^ Aii s i( IKJ ^) ~ s n{I,q)se{J,q)si-n-i{q) for all I G E A>n , J G 
E^, i > n + £ + 2p. 



Proof. For any J G E^, write J = where J G E^-p. By Lemmas [2.1| , |2.2| , we 
have (note that p is fixed) 

s e (I,q) ta 8i- p (J,q). 

Since A p > 0, for J G T,A,e- P , we can find if G Ea, p such that / = KJ G Ea,i- 
Hence when we take the sum of 7 G E^,£ on the left side of the expression, it is ~ 
to the right side summing over all J G E^_ p . This implies (i). 

To prove (ii), we fix I G E^n and take J G E^/ such that IJ G H^n+t- By 
Lemmas |2.2|, |2.1|, we have 



s n+e (IJ,q) w s n (I,q)s e (J,q). 

Thus 

^s n+ ^(/J,g) s n (I,q)s e (q). 
J 

For the reverse inequality we note for any J' G E^-p, there is if G E^ iP such that 
lif J 7 G S^n+f and 

s n+e (IKJ') ps s n (I,q)s p (K,q)st- p {J',q) ~ s n (I,q)s^ p (J',q). 
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Therefore summing over the above J 7 , we have 

2j s n+e (IJ, q) ^ 2J Sn+i{IKJ\ q) ps s„(I, q)s l - p {q) « s n (I, g)s<?(g) 

(we make used of A p > as in (i)). This proves one of the ~ in (ii). The remaining 
part follows from the same argument. 
To prove (iii), we first observe that 

2J Si(IKJ,q) w 2J s n (I,q)s i - n -t(K,q)s i (J,q) 

K: IKJe^ A ,i K: IKJeT, A>i 

4 s n (I,q)si(J,q)Si^ n ^(q) 

On the other hand, for any i^i G HA,i-n-t-2pi there exists K 2 , K 3 G XU,p such that 
IK2K1K3J G Therefore 

^2 Si(IKJ,q) ^ s n (I,q)s e (J,q) s {Kx,q) 

« s n (I, q)st(J, q)si- n -g-2p(q) 
« s n (I,q)se(J,q)si- n -e(q). 



□ 



Lemma 2.4. For a /zxed gel, 
(i) s^ +n (g) « s*(g)s„(g). 

(ii) s n (g) ~ exp(nP(g)) where P(q) is the pressure function defined in (1.1). 



Proof. From Lemma |2.3| (ii), there exist C, C > such that 

C s t (q)s n (q) < s e+n (q) < Cs e (q)s n (q), 

which proves (i). To prove (ii), we can write Csg +n (q) < (Csi(q)) (Cs n (q)). Hence 
the subadditivity property implies 

P(q) = hm lQg( ^ (9)) = inf lQg(Cg " (9)) , 

n— «x> n n n 

so that C -1 exp(nP(g)) < s n (q). The reverse inequality follows from a similar argu- 
ment . □ 

For each integer n > 0, let B n be the cx-algebra generated by the cylinders [/], 
/ G S,4,n- We define a sequence of probability measures {v n , q } on £>„ by 

^([/]) = ^f V/GS A „. (2.3) 
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Then there is a subsequence {vn k ,q}k>i converges in the weak-star topology to a 
probability measure v q . The following assertion shows that v q has the "Gibbs prop- 
erty" . 



Lemma 2.5. For a fixed q G R, ^([-Tl) ~ s n{I, q) exp(—nP(q)) for all n > 0, I G 
Proof. Let p be such that A p > 0. For any / G n and £ > n + p, we have 

■*.([/]) = E ■%.<["]> = E ££ ^r 

~ s n (J, g) ^ ? 5 (by Lemma |2.3| (ii)) 

~ s n (7, g) exp(— nP(q)). (by Lemma |2.4|) . 

Letting £ = n fc j oo, we obtain the desired result. □ 



Proof of Theorem Fix gel. Let \x q be a limit point of the subsequence 

of {i(v, + i/,o a -1 + . . . + v q o cr"^" 1 ))} in the weak-star topology. Then [i q is a 
cr-invariant measure on E^. We have for each / G n and £ > p, 

^oa-\[I\) = Yl ^([ J/ ]) 

~ s n+ e(JI, q) exp(— (n + £)P(q)) ( by Lemma 2.5) 

~ si(q)s n (I, q) exp(— (£ + n)P(q)) (by Lemma [2.3| (ii)) 

« s n (I, q) exp(— nP(q)) (by Lemma p^) . (2.4) 

This proves that p 9 is a Gibbs measure. In what follows we prove that [i q is ergodic. 
First we show that there is a constant C > such that for each I G E„, J G E^, 

ft, fcE^([ J ] n ^([ J ])) >C N ([I])M)- (2-5) 



Since \i q is supported on S^, it suffices to prove (|2.5|) for / G S^^ and J G S^. 
Note that when i > n + 2p, 

na-([J])) 

if: IKJeS A ^ +e 

^ £ s i+e (IKJ,q)exp(-(i + £)P(q)) 

K: IKJe^A,i+e 

« s n (I,q)s e (J, q)Si- n (q) exp(-(i + £)P(q)) (by Lemma ^3| (hi)) 

« s n (/,g)s<?(J,g)exp(-(n+l)P(g)) 

« ^([/]K([j]) 

from which ( j2.5| ) follows. Since the collection {[/] : / G £ n , n G N} is a semi-algebra 
that generates the Borel cx-algebra on E, a standard argument (e.g., see the proof 
of [W, Theorem 1.17]) shows that for any Borel sets A,BcS, 

- fc-i 

lim -^^(An^(B)) > Cfi q (A)fi q (B). 

i=0 

This implies that for any Borel sets j4,B C H with fi q (A) > 0, fi q (B) > 0, there 
exists n > with fl cr~ n (P)) > 0. By [W, Theorem 1.5], \i q is ergodic. 

For the uniqueness we recall that any two distinct ergodic measures must be 
singular to each other; but the Gibbs property (1.2) implies that any two fi q must 
be absolutely continuous to each other. Hence \i q must be unique. □ 

Corollary 2.6. Let \x q be the Gibbs measure in Theorem 1.1. There exists C > 
such that for any I G XU,n, J £ ^A,£ I J G T, Atn+ i, 

C-^ q ([IJ)) < » q ([I))» q ([J)) < C N {[I]) N {[J}). 

Proof. We have seen from the proof of Lemma 2.3 that for the above /, J, 
s n+i (IJ,q) « s n (I,q)s e (J,q) and from Lemma 2.4, s i+n (q) « s n (q)sg(q). By the 
definition of z/„ )? , we have 

Pn,q{[U]) ~ Vn,q([l\)Vn*([J]) 

which implies that the Gibbs measure \x q has the same property. □ 

The above property is called quasi- Bernoulli property by Heurteaux [H] (we re- 
mark that Heurteaux only introduced and studied it for measures in the full shift 
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space S). To prove Theorem |L~2"| , we need a result in [H]. Let r\ be a probability 
measure on S. For gGl, let T v (q) be the L^-spectrum of n, i.e., 

n-^oo log m _n 
where the summation is taken over all I G S n with 77 ([i]) 7^ 0. 

Proposition 2.7. ([||, Theorem 2.1]) Lei rj be a probability measure on E. Assume 
that there exists a constant C > siic/i t/iat 

rK[/J]) < Cr ? ([/])r ? ([J]) V / G £„, J G E< . (2.6) 

log 77 (/ n (a;)) 

T/ien r'(l) exists if n is a Young measure (i.e., lim — - — = constant for n 

1 ' ' n-*co logm n 

almost all x = (j,) G £, here I n (x) = [j\ . . . j n ]). 



Proof of Theorem |1.2| . For each q, let \i q be the corresponding Gibbs measure 
in Theorem 1.1 . We can view fi q as a measure on S. For t G R, let t m (£) be the L*- 
spectrum of \i q . Since /i g has the Gibbs property, it is easy to show by the definition 
of L*-spectrum that 

tP( q) - P (qt) 
logm 



W*) = TZZZZ • I 2 - 7 ) 



Note that /i 9 satisfies the condition ( |2.6|) . Since /x g is ergodic on S, it is a Young mea- 

— log fi q (l n (x)) 

sure by the Shannon-McMillan-Brieman theorem (i.e., lim equals 

n^oo n 

the entropy of fi q (with respect to a) for //^-almost all x = (ji) G £ and 7„(x) = 

b'i • • - in])- Hence by Proposition 2.7, r^(i) is differentiable at t = 1. This implies 
that -P(<?) is differential at any fixed q 7^ 0, and 

p ,, v _ P(g) -logm-r^(l) 

□ 



Proof of Theorem IAS. Let a = r'(q) with q 0. Let /i ? be the corresponding 
Gibbs measures in Theorem [Ll|, then (2.7) implies that 

qP{\) - P{q) 



E(a) = ( x G X : lim 



logm 

and 

log/ii([xi • • -Xn]) P(l)-a 



" — x logm _n logm 
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By [BMP, Theorem 1] or [LN2, Theorem 4.1], we have 

dim H E{a) < ( 1>(l) '" )<! r n (q) = ^-j— ( -aq ■ P(q)) V q 6 R. 



logm 



logm 



(2.8) 



For the reverse inequality, we see from the proof of Theorem |T72] that (1) exists 
and 

r , (1) = -gi"(g) + P(g) 
^ > logm 

By [N], we have for \i q almost all I — {%\ . . .i n . . .) G E, 

log/i,([ii, . . . ,*„]) , , 1N -qP'{q) + P{q) 



lim 

n->oo logm n 

This implies that 



lim 

n— >oo 



log ||M(x)M(ax) ■ ■ • M^^x 



71 



Therefore we have 



dim// £"(«) > dim// 



logm 

P\q) = a ji q — a.a. x G S. 

-ga + P(g) 
logm 



□ 



3. Nonnegative matrices 

In this section, we always assume that M is a function on taking values in the 
set of all d x d non-negative matrices and satisfies (HI) and (H2) defined in Section 
1. Let q > be fixed. Then s n (I,q) and s n (q) in (2.1) are reduced to 

s n (I,q) = ||M/|| 9 V/GS Ajn and s n (q) = £ s n (/,g). 



For convenience, we let 



mm mm 

s,te{l,2,...m} l<i,j<d 



E E M H • ( 3 -^ 

fc=i i^es A , fc „, t 



Then b > by (H2). 

We will reformulate the three theorems in the previous section. The proofs are 
almost the same and for simplicity, we only point out the differences. Here Lemmas 



2~T| , |2.2| do not hold anymore; on the other hand we can use (H2) to replace these 



lemmas to obtain an analog of Lemma 2.3 



li 



Lemma 3.1. For a fixed q > 0, 

(i) s e+1 (q) « SfXq). 

(ii) For I G £ A ,n, J2jS n+e (IJ,q) « ~ s n (J,g)s*(g) w/jere toe 
/irs£ (second) sum is taken over all J G XU,£ s,uc ^ £ SA, n +£ ( JI E T,A, n +e 
respectively). 



(iii) EfcLi Ex: /xj e s A , l+fe s i+k {IKJ,q) « s n (J, g)s*( J, q)si- n - t (q) for all I G 
S Aj „, J G i > n + & 

Proof. For any / G Ea,m-i, write I = iJ with J G XU,£- Using ||M/|| < ||Mj||||Mj||, 
we have for q > 0, 

s e+1 (q)<m( sup ||M;|| 9 )s £ (g). 

i£{l,2,...,m} 

That is, Sf + i(q) =4 S K<?)- F° r the reverse inequality, since for any J G E^, 

r r 

E E 11^11 = IKE E m^mjH^&iimjH, 

fc=i k&ax- KJe^ A ,e+k k=i xeE Ajfc : KJez A , e+k 

it follows that 

E E ii^i* > (y^nw- 

This combines with S£+i(q) =4 sg(q) imply that 

r 

s t+i(o) > ^2s e+k (q) > s e (q) 
k=i 

and completes the proof of (i). 

To prove (ii), it follows from \\Mu\\ < ||M/||||M/|| that Y.j& AX .u& A>n+l s n+e(H,q) 
=4 s n (I,q)se(q). For the reverse inequality, we use (H2) as above to conclude that 
for any J G T, A/ , 

r 

E E ii m ^ii > ii M i ii ii m j ii • 

k=l K£Y,A,k'-IKj£Y,A,n+e+k 

Hence ELi E^ 6 E A , fc :7^j e s A , n+ , + , II^jII 9 ^ II M i II 9 II M J II 9 and therefore, summing 
over the J G Sa/, 

r 

E E s„+£+fc(/£,g) ^ s n (I,q)s e (q). 

k=l L^Y. A,e+k ]I ie E A,n+e+k 

Since 

E s n+e (IJ,q)^ E s n+e+1 (IJ,q), 

j£^A/-U'^<A,n+l ^E/1,(+I : ^Ea,ii+<+1 
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we have 

r 

(IJ,q) fc= s n (I,q)st(q). 

A,n + l A,n+l+k 

This completes the proof of an rs in (ii); the other « follows from an identical 
argument. 

To prove (iii), we have, 

Y Y s i+k (IKJ,q) < s n (I,q)si(J,q)Y s i+k-n-e(q) 

k=l K:IKJ€T, A ,i+k k=l 

~ s n (g)s^(J,g)si_ n _ £ (g) (by {%)). 

On the other hand, for any W G T,A,i- n -£, by (H2), there exist 1 < &i < r, K\ G 
Syi^j such that IKiW G X^j^+j^ and 

„ MI-HI II 

W > " r "" 1' , 

where 6 is defined by (|3.1|). By using (H2) again, there exist 1 < k 2 < r, K 2 G S^ )fc2 



such that IKiWK 2 J G Si ii+ i 1+ fc 2 and 



| /A - lM --AVII>^™> 62||/|ll|W ' lll|J| 



ELi mfc (ELi mfc ) 2 

Therefore we have 

Y Y s i+k (IKJ,q) )p s n (I,q)st(J,q)Si- n -i(q) 

fc=l K:IK.J&Y, A , l+k 

and (iii) follows. □ 
We now state the corresponding theorems as in Section 1. 

Theorem 3.2. Suppose M is a function on T,a taking values in the set of all d x d 
non-negative matrices and satisfies (HI) and (H2). Then for any q > 0, there is a 
unique Gibbs measure \i q on Xa as in Theorem 



The proof is almost identical with that of Theorem using Lemma 3.1. The 
only adjustment is to replace 

fi q ([/] n a-*([J])) « n q (I)n q (J), i>n + 2p. 

by 



2r 



J> 3 (M n <r*([J])) « ^(/)/i 3 (J) vo«. 
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We use the same proof as in Section 2 for the next two theorems. 



Theorem 3.3. Under the same conditions of Theorem \3.3j, P(q) is differentiable 
for any q > 



Theorem 3.4. Under the same conditions of Theorem \3. \ we have for any a = 
P'{q), q > 0, 

dim H E(a) = (—ocq + P(q)) 

logm 

where E(a) = { J = (ji) G S A : lim™-^ log \\M h ■ ■ ■ M jn \\/n = a}. 

To relate Theorem 3.4 to the classical random product of matrices, we let {Y n } be 
the i.i.d. random variables that take values M\, . . . ,M m , invertible matrices and 
with uniform distribution, then lim^oo - log \\Y n . . .Yi\\ = A a.s. and A is called the 
upper Lyapunov exponent ([FK], [BL, Chapter 1]). In comparison with Theorem 
3.4, we let S A = E be the space of full shift (i.e., all the entries of A are 1), then 
-P(O) = logm. The limit of the random variables corresponds to the case for q = 0, 
A = -P'(O) and dim H E(X) = -P(O)/ logm = 1 (the existence of the derivative follows 
from some additional assumptions on the Mj ([BL, p. 119]). 

We remark that if S A = S, then condition (H2) is reduced to a more simple form: 
YTk=i H k > where H = Mi + - • • + M m . The condition is essential for the theorems 
in Section 3. Indeed we have 

2 \ / 2 

Example 3.5. Let M x = ( n ] and M 2 = ( ) . Then // .\/, • 

is reducible. Since Mj — ( | where \J\ = n and k is the number of 2's 

/|J|=n II 

^^W + 3 fc ) 9 >max{^^W ^ tf\ 3 fcg } = max{2^ +1 ), (1 + 3*)"} 

k=0 ^ ' k=0 ^ ' k=0 ^ ' 

and 

(?) (2 n + 3 h ) q < (fc) 2 ^ 2 "* + ^ = 2<? (2 n( " +1) + (! + 3 T)- 

fc=0 ^ ' fc=0 ^ ' 

We have = max{(g + 1) log 2, log(l + 3 9 )}, which is not differentiable at q — 1. 
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2 







:) 


2" 








3 fe 




l~ik=< 



We see that the Gibbs measure p q in Section 2 has the quasi-Bernoulli property. 
However for the case of non-negative matrices, only p q ([IJ]) < C p q ([I])p q ([J]), 
I E £„, J G Hg. The following example shows that the reverse inequality may not 
hold. 

Example 3.6. Let M 1 = ^ q ^ ^ anc ^ ^2 an arbitrarily positive matrix, then 

H — Mi + M 2 is an irreducible positive matrix. Let J = 1 ... 1 (n-times), then 
\\Mj\\ = n + 2 and hence 

||Mj|| \\Mj\\ > -\\Mjj\\ . 

Since p q {I) ~ exp(— nP(q)) ■ ||Mj||, / e S n , we see that there does not exist C > 
such that C'n q ([J])n([J]) < fi([JJ]). 



4. Examples and remarks 

Consider the classical Bernoulli convolution X = (1 — p) _1 X^iP n ^« where the 
X n 's are i.i.d. random variables which take values 0,1 and with probability 1/2 
on each value. Let \i p be the distribution measure of X. It is well known that for 
< p < 1/2, the measure is a Cantor type measures. It was proved recently that fj, p 
is absolutely continuous for almost all 1/2 < p < 1 [S], however, it is still not clear 
which p p is absolutely continuous or singular. The question has been subjected to 
intensive investigation, the reader can refer to the survey articles [L], [PSS] and the 
references there. The interest of the Bernoulli convolution in our setting is that the 
p p satisfies the self- similar identity 

1 1 

ftp — ^p^i + 2^ p ^ 2 

where Six = px, S 2 x = px + (1 — p); {Si, S 2 } is the iterated function system (IFS). 
The support of p q is [0, 1]. For < p < 1/2, the Sj(0, l)'s are disjoint (as in the 
basic intervals of the Cantor set); for 1/2 < p < 1, the Sj(0, 1) overlaps which is the 
source of difficulty 

For p = (\/5 — l)/2, the reciprocal of the golden ratio, it was shown by Erdos 
that p p is singular. In order to consider the multifractal structure of p p , we can put 
the IFS {Si, S 2 } to a new set of IFS {i?j}f =1 with no overlap: 

Ri(x) = p 2 x, R 2 (x) = p 3 x + p 2 , and Rs(x) = p 2 x + p. 

Then the measure p p satisfies 

p p {[ii...i n \) w ^\\M h ...M in \\ 
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where [h . . . i n ] = Rh ■ ■ ■ #i„([0, 1]) and 

Ml = (°0' M2= {\\) a " d Ma= 0°) 

([LN1], [Fel], [FeO]). The {M;}f =1 satisfies the conditions (H2). Hence by Theorem 
3.2, 3.3, P(q) is differentiable for q > and the multifractal formalism holds. 

Actually more can be said about the L 9 -spectrum r(g) of fi q : an explicit formula 
was given in [LN1] for q > and was extended to q < in [Fel]. By using the 
formula it was proved that r(q) is differentiable (actually real analytic) on K. except 
one point in M~ in [Fel]. 

The above example of Bernoulli convolution gives rise to another interesting ques- 
tion. Note that the above example is a special case of the overlapping IFS that can 
be reduced to new sets of IFS with no overlap and the calculation of the r(q) can 
be converted into the product of matrices as in (1.4). Such IFS forms an important 
subclass of those that satisfy the weak separation condition ([LN2], [LNR])(it will 
be interesting to classify this subclass of IFS). Under the weak separation condition 
it was proved that the multifractal formalism is valid provided that r(q) is differen- 
tiable [LN2]. However we do not know its differentiability in the general case. In a 
forthcoming paper [Fe2], the first author proves the differentiability of r(q) for q > 
in the case that the IFS is equicontractive and satisfies the finite type condition (see 
[NW]). 

The behavior for q < is also important for the multifractal analysis. There is 
no problem when M is a positive matrix-valued as we consider in Section 2. For 
the non-negative matrix-valued M, Mj can be 0, we have to modify the pressure 
function P(q) in (1.6) slightly: 

P(q) = lim -log V \\Mj\\ q (4.1) 

JeAf n 

where M n consists of all the J G XU,n such that Mj ^ 0. It is clear that if Mj ^ 
for all J G £A,n) then the super-additivity of the sum in (4.1) implies that the limit 
exists. We include a simple proposition with = E to set up the consideration: 

Proposition 4.4 Suppose M 1: • • • , M m are non- negative matrices and H = YlT=i 
is irreducible, then the limit in (4.1) exists for each q < . 

Proof. By the irreducibility, there exists integer r with J2k=i H h > 0. Hence 
there is a constant C > such that for any two finite indices /, J, there exists 
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K o e ULi s fc satisfying 

0< \\M IKoJ \\ KCWMjWWMjI (4.2) 

Denote by s n = Y..jeN n \\ M j\\ q - Then ( 4 -2) implies s n s £ < C~ q YJk=i s n+t+k- From 
(4.2) we also deduce that for any finite index /, there exists i G Si such that Mn ^ 0; 
Since ||M/j|| < Ci||M/|| for some constant C\ > 0, we have s n < C^ q s n+ i for any 
integer n, £. It follows that s n st < C's n+ i +r for some constant C > (depending on 
q), which implies that a n = ^-Sn-r is super-multiplicative. This yields the existence 
of the limit. □. 

The differentiability of such P(q) for q < is still unknown. We know that in 
the above Bernoulli convolution of the golden ratio, it is possible for the P(q) to be 
non-different iable at a point of q < [Fel]. On the other hand, it is known that 
by imposing some stronger conditions on the matrices, the pressure function P(q) 
is analytic near q = (see e.g., [BL, Theorem 4.3]). 

Finally we remark that we do not know whether the theorems can be extended 
matrices with entries in R. An important theorem concerning this is in [BL, Theorem 
4.3] for the analyticity of P(q) near zero. Much closer to our development is the 
scaling functions: f(x) = Y^o c if(^ x ~ ^ ^ s known that such function can be 
put into matrix form as in the previous example [DL1]. Daubechies and Lagarias 
studied the multifractal formalism of the well known scaling function _D 4 [DL2]. 
They showed the differentiability of the corresponding r(q), but the consideration 
depends on the two 2x2 matrices involved to have a common eigenvector. There 
are some extensions in [LM]. 
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